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Abstract: We are interested in the modehng of wave propagation in poroelastic media. We 
consider the biphasic Biot's model in an infinite bilayered medium, with a plane interface. 
We adopt the Cagniard-De Hoop's technique. This report is devoted to the calculation of 
analytical solutions in two dimensions. The solutions we present here have been used to 
validate numerical codes. 
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Solution analytique pour la propagation d'ondes en milieu 
poroelastique stratifie. Partie I : en dimension 2 

Resume : Nous nous interessons a la modelisation de la propagation d'ondes dans les milieux 
infinis bicouches poroelastiques. Nous considerons ici le modele bi-phasique de Biot. Cette 
partie est consacree au calcul de la solution analytique en dimension deux a I'aide de la 
technique de Cagniard-De Hoop. Les solutions que nous presentons ici ont ete utilises pour 
la validation de codes numeriques. 

Mots-cles : Modele de Biot, ondes poroelastiques, solution analytique, technique de Ca- 
gniard de Hoop. 
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Introduction 

Many seismic materials cannot only be considered as solid materials. They are often porous 
media, i.e. media made of a solid fully saturated with a fluid: there are solid media perforated 
by a multitude of small holes (called pores) filled with a fluid. It is in particular often the 
case of the oil reservoirs. It is clear that the analysis of results by seismic methods of the 
exploration of such media must take to account the fact that a wave being propagated in such 
a medium meets a succession of phases solid and fluid: we speak about poroelastic media, 
and the more commonly used model is the Blot's model [H El E]. 

When the wavelength is large in comparison with the size of the pores, rather than regarding 
such a medium as an heterogeneous medium, it is legitimate to use, at least locally, the theory 
of homogeneization [UlIS]. This leads to the Blot's model [HEIE] which involves as unknown 
not only the displacement field in the solid but also the displacement field in the fluid. The 
principal characteristic of this model is that in addition to the classical P and S waves in a 
solid one observes a P "slow" wave, which we could also call a "fluid" wave: the denomination 
"slow wave" refers to the fact that in practical applications, it is slower (and probably much 
slower) than the other two waves. 

The computation of analytical solutions for wave propagation in poroelastic media is of high 
importance for the validation of numerical computational codes or for a better understanding 
of the reflexion/transmission properties of the media. Cagniard-de Hoop method [Sj [7] is a 
useful tool to obtain such solutions and permits to compute each type of waves (P wave, S 
wave, head wave...) independently. Although it was originally dedicated to the solution to 
elastodynamic wave propagation, it can be applied to any transient wave propagation prob- 
lem in stratified medium. However, as far as we know, few works have been dedicated to the 
application of this method to poroelastic medium. In [11] the analytical solution to poroelas- 
tic wave propagation in an homogeneous 2D medium is provided. 

In order to validate computational codes of wave propagation in poroelastic media, we have 
implemented the codes Gar6more 2D [9j and GarGmore 3D [10] which provide the complete 
solution (reflected and transmitted waves) of the propagation of wave in stratified 2D or 
3D media composed of acoustic/acoustic, acoustic/elastic, acoustic/poroelastic or poroelas- 
tic/poroelastic layers. The codes are freely downloadable at 



|http : //www . spice-rtn . org/library/sof tware/Gar6more2D 

and 



http : //www . spice-rtn . org/library/sof tware/Gar6more3D[ 



We will focus in this paper on the 2D poroelastic case, the 2D acoustic/poroelastic case 
is detailed in [8] and the three dimensional cases will be the object of forthcoming papers. 
The outline of the paper is as follows: we first present the model problem we want to solve 
and derive the Green problem from it (section 1). Then we present the analytical solution to 
the wave propagation problem in a two-layered 2D poroelastic (section 2). Finally we show 
how the analytical solution can be used to validate a numerical code (section 3). 
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1 The model problem 

We consider an infinite two dimensional medium (Q = R ) composed of two liomogeneous 
poroelastic layers ri"*" = Rx] — oo,0] and fi" = R x [0, +oo[ separated by an horizontal 
interface T (see Fig. [1]). We first describe the equations in the two layers ( §l.ip and the 
transmission conditions on the interface F ( §1.2|) . then we present the Green problem from 
which we compute the analytical solution ( §1.30 . 



First Layer 
y = ' 



Second Layer 



Figure 1: Configuration of the study 



1.1 Poroelastic equations 

We consider the second-order formulation of the poroelastic equations [H EJ [3] : 

in nx]0,T], 

in nx]0,T], 

in nx]0,T], 

P + f3V -Us + V-W = F„, mnx]0,T], 

m 

Usix, 0) = 0, W{x, 0) = 0, in Q, 

tlsix, 0) = 0, W{x, 0) = 0, in Q, 



Pftjs + p^nW + ^W + VP = F^., 

S = AV ■ [/, /2 + 2p£{Us) -f3Pl2, 
1 



with 



(1) 



dT,- ■ 
(V • S)j = y^ ''' y i = 1,2. As usual I2 is the identity matrix of ^^2(2^)) 



and eiU s) is the solid strain tensor defined by: 
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In ([T]), the unknowns are: 

• Us the displacement field of solid particles; 



W = 4>{U f — Us), the relative displacement, Uj being the displacement field of fluid 
particles and (p the porosity; 



• P, the fluid pressure; 

• S, the solid stress tensor. 

The parameters describing the physical properties of the medium are given by: 

• p = (j) pf + {1 — (t>)ps is the overall density of the saturated medium, with ps the density 
of the solid and pj the density of the fluid; 

• Pw = o-pf/cj), where a is the tortuosity of the solid matrix; 

• /C = k/i], where k. is the permeability of the solid matrix and r/ is the viscosity of the 
fluid; 



m and P are positive physical coefficients: P = 1 — K^/Ks 

and m = [(p/Kf + (/3 — (l))/Ks]~ , where Ks is the bulk modulus of the solid, Kf is the 

bulk modulus of the fluid and Ki, is the frame bulk modulus; 



• /x is the frame shear modulus, and X = Kb — 2p/3 is the Lame constant. 

• Fu, Fw and Fp are the force densities. 

To simplify this study, we consider only the case of a compression source 

Fu{x, y, t) = fuV{6r, 5y-h) f{t) and F^(x, y, t) = U^{6^ 6y-h) f{t) 

and a pressure source Fp = fp6x Sy~h fit), where /«, f^ and fp are constant and / is a regular 
function source in time. We can generalize this approach for other types of punctual sources 
such as for instance 

Fu = fu^ X {6x5y-h) fit) and F^ = f^jV x {d^^Sy-h) f{t). 
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1.2 Transmission conditions 

Let n be the unitary normal vector of T outwardly directed to 0,^ . The transmission condi- 
tions on the interface T between the two poroelastic medium are [6j: 

r ut = u-, 

W+ ■ n = W ■ n, 
P+ = P-, 



(2) 



1.3 The Green problem 

We won't compute directly the solution to ([2]) but the solution to the following Green problem: 



^±--± , „± ".i 



1 



Pj Uf + Pw W + p^ W + "^P = fvN{5x Sy-h)St 

a^ = A^V ■ufl2 + 2p^e{uf) - (3^ p^ I2 
J_p± + /3± V . w± + V • «;± = fp6x 6y^h f{t) 



inQ=^x]0,r], 

inS7=^x]0,T], 
in[7=^x]0,r], 
inO=^x]0,r], 



m7=u+, onrx]0,r] 

w n = w^ n, onTxJOjT] 

p- =p+, onrx]o,r] 

a'n = a+n, onrx]0,r]. 



(3a) 

(3b) 
(3c) 
(3d) 

(3e) 
(3f) 

(3g) 
(3h) 



The solution to ([T]) is then computed from the solution to the Green Problem thanks to a 
convolution by the source function. For instance we have: 

P^ix,y,t) =p^{x,y,.)* f{.) = / p+{x,y,T)f{t-T)dT 

Jo 

(we have similar relations for the other unknowns). We also suppose that the poroelastic 
medium is non dissipative, i.e the viscosity rj^ = 0. Using the equations (|3c|3d|) we can 
eliminate a and p in ([3]) and we obtain the equivalent system: 

' p^iif + pf iv^ - a^ V(V ■uf) + p^V X {V X uf) - m^l3^V{V ■ w^) 

= ifu - P^m+fp)V{6x 6y.h) 6t, (4) 

t pfut + pt w^ - m^P^ V(V • uf) - m± V(V ■ w^) = {U - m+fp)V{6x 6y^h) St, 

with a~ = A^ + 2p^ + m~/3~ . 
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And the transmission conditions on T are rewritten as: 



'^sy ~ '^sy^ 



W, 



W' 



y y ' 

m+/3+ V • W^ + 771+ V • w^ = m~ P~ V • u^ + m^ V • w 



(5) 



(A- + m+/3+ )V • M+ + 2/i+9j^n+ + m+/5+V 



w 



(A + m /3 )V • li^ + 2/i 9j/Us« + "^ /^ V • it» . 



We split the displacement fields u^ and w on irrotational and isovolumic fields (P-wave and 
S-wave) : 



u 



± 



VQ- + V X *= 



w 



^ = vet + vx^t. 



(6) 



We can then rewrite system Q in the following form: 

' A+Q+ - B+Ae+ = SJy^h^tF, mn+x]0,T] 
A-O' - B-AQ- = 0, mn-x]0,T] 



(7) 






4^^ 

Pw 



where G± = (9^, 6^)*, F = (/„ - /3+m+/p, /^ 
matrices: 



inS7^x]0,r] 
1^^ fpY^ ^^ cind S^ are 2x2 symmetrical 



A 



± 



P/ P^ 



± 



and 



B 



Vs^ 



A±+2/i±+m±(/3±)2 m±/3± 



±/o± 



m P 



m 



± 



IJipt 



\ P^pw - Pf 



±2 



is the S-wave velocity. 



We multiply the first (resp. the second) equation of system ([7]) by the inverse of A'^ (resp. 
A-). The matrix A+~^B+ (resp. A-'^B') is diagonalizable: A^~^B^ = V^D^V^~^, 
where V is the change-of-coordinates matrix, D = diag{Vpr , Vp^ ) is the diagonal matrix 
similar to A^ B^, Vpr and Vp^ are respectively the fast P-wave velocity and the slow P- 
wave velocity (Vp^ < Vpr). 



Using the change of variables 



$= 



i^pp^PsY 



v^-'e^, 



(8) 
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we obtain the uncoupled system on fast P-waves, slow P-waves and S-waves: 



$+ _ D+A$+ = 6rA,-hS, F+, in Q+x]0, T] 



^- - D-A<^- = 



^xOy-hOt 

0, 



^ 



± 



4*^ 

Pw 



in n-x]0,T] 
in O=^x]0,T] 

in 0±xlO,Tl 



with F+ = {A+V+)'^F = (F+^,F+J*. 

Finally, we obtain the Green problem equivalent to ^: 



^f-V+^A^+ = 5,,6y_hStF+, ie{Pf,Ps} y>0 



$+ - ^5 A$+ 







$; 



r-2 



0, 



Bi^pp'^ps^^l'^pp^ps^'^s: 



y>0 
y<0 
y = 



(9) 



(10) 



where we have set ^^ = ^^ in order to have similar notations for the Pf, Ps and S 
The operator B represents the transmission conditions on T: 



waves. 



B 






Vt.dy Vt,dy 



a 



Pf 



Pw 



1351 

Bqi 



^42 
B52 
Bq2 







P^{dyy-d': 



,2 N 
XX J 



V7,d^ -pr,a 



11 "^a;' 



-^iidy 



Bu 

B5A 

Bqa 



\2'-'x. 



-^12 C'j/ 



-^21 5y -^22 ^S/ 



^45 
^65 



-5y 
Pw 



2p-dly 





\^IA 




^h 




^^ 




^pf 




•^Ps 




.^5 . 
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where Pj ■ , i,j = 1,2 are the components of the change-of-coordinates matrix V and 



_ m+(/?+P++P+) , . _ m+(/3+P++P+) ,. 

^41 - — ^ % , O42 - —2 %' 

m-(/3-7^n+7^,l) 2 . ^ m-(r7^f2 + ^22) ^2. 

O44 — ^-2 C'ti > •^45 — 32 "^it) 

B,i = 2A.+P+ d^,y ; ^52 = 2/.+P+ d^y ; ^54 = -2^-7^11 d^y ; ^55 = -2//-Pr2 d%; 

£361 - ^:i % + 2/i l-^uOyy, 

Vpj 

(A+ + .n+/?+^)P++m+r^2+2 ,.2 , g,+^+.2. 
^64 = ^ ^ " ^^^^2' ' "^ " ''' ^'f + '^l^"Plld'yy. 



(A- + m-r^n + "i~r ^21 ^2 , o,,-^- ^2 

;^ _ (A~ + m-l3-^)V^2 + ^~/^~^2'2 a2 , o,.-n- fl2 
065 - 32 % + ^l^ ' 12'^yy 

To obtain this operator we have used the transmission conditions ([5]), the change of variables 
©-([H]) and the uncoupled system Q. 

Moreover, from the unknowns ^pt, ^pg and $^ we can determine the solid displacement 
uf and the relative displacement w^ by using the change of variables presented below. 

2 Expression of the analytical solution 

To state our results, we need the following notations and definitions: 

1. Definition of the complex square root. For q G C\IR~, we use the following 
definition of the square root g[q) = q^'"^: 

g{qf = q and ^e[g{q)] > 0. 

The branch cut of g{q) in the complex plane will thus be the half-line defined by {q G 
IR~} (see Fig. [2]). In the following, we use the abuse of notation g{q) = i\/—q for 
qeJR-. 



2. Definition of the functions Kf . For i G {Pf,Ps,S} and g € C, we define the 

functions 

\ 1/2 



K^:=K+{q)= I ^ + 



V^ 



1 

±2 
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Kg) I 



+vr 




^e{q) 



Figure 2: Definition of tlie function x 



[x 



,1/2 



3. Definition of tiie reflection and transmission coefficients. For a given q £ C, we 
denote by TZpfPf{q), TZpfPs{q), TZpfsiq), Tpfpf{q), TpfPs{q), and Tpfs{q) tfie solution 
to tlie linear system 



Aiq) 



T^PfPfiQ) 
T^PfPsiq) 

T^pfsiq) 

Tpfpf{q) 

TpfPsiq) 
'^pfsiq) 



2k+ f{q)Vi 



Pf 



iqVl, 
-Kpfiq)V^l 



V, 



Pf 

2ig^+4/g)P+ 
(A+ + m+P+^)Vti + m+p+V^^ 



V, 



2/x+P+4/g) 



Pf 



and by TZpsPf{q), TZpsPsiq), TZpssiq), TpsPf{q), TpsPs{q) and Tpss{q) the solution to 
the linear system 



A{q) 



T^psPfiq) 

TlpsPsiq) 

T^pssiq) 

'TpsPfiq) 

TpsPs{q) 

Tpssiq) 



24s(9)^A 



m 



v-^ 



^q'Pu 
-4s(Q)'Pt2 



;{(i^vt, + n 



22) 



' Ps 

2i qfi+4^{q)rt2 
(A+ + m+/3+^)P+ + m+(3+V^^ 



V, 



2fi+Vt24sil) 



Ps 
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11 



where the matrix A{q) is defined for q £ C by: 



A{q) 



-[qV+ -igP+ 



-f^Ml) 



^Pf\H)i 11 



12 



-4fiq)V+ -4siq)'P22 



-\q 



Pf 



Aii{q) 

Ai(9) 



Ai2{q) 
ArAl) 
AQ2{q) 





A^q) 



i^^ii 19^12 

-i^PfiQ)'Pu -'^Psil)^ 12 

-Kpj:{q)V2i -'^Psil)P22 

Auiq) A45{q) 

A54{q) A55{q) 



%(«) 



Pw 



-2ig^+K+(g) ^4(9) 



A(i5iq) 





A^eiq) 
-2iqfj,-Kg{q) 



with 
Aiiiq) 

Auiq) 



m^ 



V, 



+ 2 



[/3+P++P2+1]; M2iq) 



Pf 
m 



V, 



— [mi + ^2l]; M^ii) 






[p+vt^ + n 



+1 • 
22J ' 



Pf 



V, 



[/3-K2 + P22]; 



Ps 



A^iiq) = 2iq44^{q)Vti; A^q) = 2ig^+4^(g)P+ ; ^3(9) = Ai+(4 (9) + ?'); 
A-,i{q) = 2i qn-Kpf{q)V^^; ^5(9) = '^'^QlJ''Kp^{q)P{2' Ael?) = -f^'i^g (q) + q^); 



Aeiiq) 
Ae2{q) 
Aeiiq) 
Ae5{q) 



(A+ + m+/3+^)Vl^ + m+13+V. 



21 



V, 



+ 2fi+KJ,/{q)r 



u'l 



Pf 



iX+ + m+P+^)Vt2+m+f3+V^^ , „ ^ ^ . .2 



V, 



Ps 

2\ 



(A-+m-/?-')Pfi+m-/3-P2i . - - 2 



V, 



Pf 



{X-+m-p-')V{2-m-p-p^^ , „ - -2 



V, 



+ 244^{qrV+; 
-2fi-Kpj: {q)r^^; 
+ 2fi-Kp^ {q)V^2- 



Ps 



We also denote by Vmax the greatest velocity in the medium: 

Knax = max(Fj+ , V^^, V^, Vpf, Vp^, Vg). 
We can now present the expression of the solution to the Green Problem: 
Theorem 2.1. The solid displacement in the top medium is given by 



with 



u 



uj{x,y,t)= u^{x,y,T)dT, 
Jo 



Upj: + Up^ + t^pfpf + t^pfPs + '^PfS + ^PsPf + ^PsPs + ^PsS 
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and the solid displacement in the bottom medium is given by 



with 
where 



Us {x,y,t) = \ V {x,y,T)dT 
Jo 



V — l^pfpf + ^pfPs + ^PfS + ^PsPf + ^PsPs + ^PsS-" 



u'^r is the velocity of the incident Pf wave and satisfies: 



^Pf,xi^^y^t) 



''Pf,yix,y,t) 



t^Pfix,y,t) =0. 



'Pii^Pf tx 



^ii^p/ tiy-h) 



F+/ I'KT^fW^' 



'Pf 



else 



We set here r = (z^ + (y — /i)^)^'^ and to = r/Vpf denotes the arrival time of the 
incident Pf wave. 



Pf 



Up is the velocity of the incident Ps wave and satisfies: 






tx 



i^Ps,y(a;,y,i) 



'Pt2Fps t{y-h) 



ift>to 



V+' 27:r.JW~ti 



^tsi^^y^'^) = 0- 



else 



We set here r = (x^ + {y — h)'^)^''^ and to = ^/^ps denotes the arrival time of the 
incident Ps wave. 



I'pfPf is the velocity of the reflected PfPf wave (the Pf reflected wave generated by 
the Pf incident wave) and satisfies: 

yv{t)4,j{v{t))npfPf{v{t)) 



^^pfpf.xi^^y^^) 



'PfPf,x\ 



'PfPf,y 



{x,y,t) 



4/{vit))npfPf{vit)) 



.-p+ p+ 
' 11-^ P/' 



X ^Pf 

if tf^ < t <to and — > 



r Kn 



vr 



aAF^ 



' 11-^ p/' 
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Re 



"pfpfj^^y^t) 



7(t)4f(7(i))^P/P/(7(i)) 



Re 



^PfPf^y^^^y^^) 



4/(7(t))^P/P/(7(i)) 



■niFpp 



ift>to 



' 11-^ P/' 



UpJx, y,t) = 0. e/se 

W^e set /lere r = (x^ + (y + /i)^)^'^ and to = '^/^^f denotes the arrival time of the 
reflected PfPf volume wave, 



th = iy + ^)j7^-^ + ^ 

V Pf '^inax •'max 

denotes the arrival time of the reflected PfPf head wave and the complex functions 
V := v{t) and 7 := 7(t) are defined by 



vit) = -i 



y + h 



\ 



1 t^ 2;t 

^ H — 5-1 for tfi < t <tQ and x < 0, 



Vpf 



vit) = i 



y + /i 



\ 



1 t2 3,^ 






/or tfi < t < to and x > 0, 



and 



, . . X V + h 
7(t) = -i-t + ^ 



'\ 



t' 1 



r^ y 



/or t > to- 



Pf 



I'pfPg is the velocity of the reflected PfPs wave and satisfies: 



•p+ 77+ 

^pfps,xi^^y^t) = — ^^^ 



^P/Ps,j/(^'y'*) = ^^^ 



yv{t)npfps{v{t)) — {t) 



du 

K+^{v{t))npfPs{v{t)) — {t) 



iftfi<t<to 

and \Qm (t;(to))| > 



Kr 



''pfPsJ^'y'^) 



vt.F^ 



+ 
i^"^ Pf 

vr 



Re 



+ ^ j.\ ^"^ PffO 

^PfPs^yi^^y^t) = ^e 



yi{t)npfPs{i{t))^{t) 



K+^{j{t))npfPs{-f{t))^{t) 



ift>to 



^'pfPs^^^y^^) = 0- 



else 



Here to denotes the arrival time of the reflected PfPs wave (its calculation is similar 
to the calculation of the arrival time of the transmitted wave, see the appendix of ^8]) 
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and tf^ denotes the arrival time of the PfPs head wave: 



h = y\ 



y^s ^^n 



' ^kU 



1 \x\ 

+ 



V p r max "^ max 



For tf^< t < to, the function v{t) is implicitly defined as the only root of 



1 



1/2 



1/2 



qeC^J^{q,t)=y\^ + q'\ +h 



V, 



Ps 



Vr 



+ 



+ iqx — t, 



Pf 



such that Qm f x ^^^ ' 1 < 0. 

For t > to, the function 7(t) is defined as the only root of q £ C >-^ J^{q,t) whose real 
part is positive. 

^PfS ^^ ^^^ velocity of the reflected PfS wave and satisfies: 



^PfS,xi^^y^^) 






3f?e 



vr 



4{v{t))npfs{v{t))—{t) 






di) 

vv{t)npfs{v{t))—{t) 



iftf^<t<to 

and |9m (f (to))| > 



V^ 



^pfsJ^^y^*) 



^pf 



die 



vr 



4iiitmpfs{iit))^it) 



^PfS^yi^^y^'^) 



F+ 



^e 



vr 



i7(t)7^P/5(7(i))^(i) 



ift>to 



Vpjsix,y,t) = 0. 



else 



Here to denotes the arrival time of the reflected PfS wave and t^ denotes the arrival 
time of the reflected PfS head wave: 



h = y\ 



1 



T/+2 y2 

V c m 



1 +/> ' ' 



1 \x\ 

+ 



F+/ K?.ax Vr 



For tf^<t< to, the function v{t) is implicitly defined as the only root of 



1 



1/2 



1 



1/2 



qeC^J^{q,t)=y\^ + q'] + h I ^-^ + q^ ] + iqx - t, 



Vo 



V, 



Pf 



such that 9m ( x "^^ ' I < 0. 

For t > to, the function ^{t) is defined as the only root of q €z C ^^ J-{q,t) whose real 
part is positive. 
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^~PsPf ^-^ ^^^ velocity of the reflected PsPf wave and satisfies: 



^PsPfA^^y^i) 



PsPf,y 



TT 



'x,y,t) 



' n^Ps 



Re 



TT 



dx) 

■iv{t)npsPf{v{t))—{t) 



K+pAv{t))npsPf{v{t)) — {t) 



dt 



iftj^<t<to 

and 1 9m (t;(to))| > 



V„ 



^psPf,x{^^y^t) 
^^psPfJ^^y^t) 



IT 



-^^^^Re 

TT 



n{t)npsPf{-f{t))^{t) 



4f(7(i))^P.P/(7(i))Sw 



dt 



ift>to 



^tspA^iVi't) = 0- 



else 



Here to denotes the arrival time of the reflected PsPf wave and t^ denotes the arrival 
time of the reflected PsPf head wave: 



1 1,1 1 \x\ 

tfi = yj — - T77^ \-h, I — - T7?. \- 



y+/ K^ax 



V r,„ "^ max ^ max 



Ps 



For ti^ < t < to, the function v{t) is implicitly defined as the only root of 



1 



1/2 



1 



1/2 



qeC^J^{q,t)=y\^-^ + q'] + h I —-^ + q' ] + iqx - t, 



V, 



Pf 



T/+2 

^Ps 



such that 9m I x ^^^ I < 0. 

For t > tQ, the function ^{t) is defined as the only root of q £ C i-^ J^{q,t) whose real 
part is positive. 

• ^~PsPs ^-^ *^^ velocity of the reflected PsPs wave and satisfies: 



''PsPs,xix^y^i) 
^PsPs,yi^^ y. 

^psPsJ^^y^t) 

^PsPs,y(.X^ y. i) 
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iv{t)K'^f{v{t))TZpsPs{v{t)) 



T^Vtf^ 



,-p+ p+ 



KJ,f{v{t))npsPs{v{t)) 



X V 
ifth<t<to and - > ^"^ 



r V„ 






vr 



aAF^ 






Re 



4>s {l{t))'RpsPs{l{t)) 



TT.JW—tl 



.-p+ TP+ 



ift>to 



else 
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We set here r = {x^ + {y + h)'^)^''^ and to = ^/^ps denotes the arrival time of the 
reflected PsPs volume wave, 



th = {y + h) 



1 \x\ 

+ 



Vpr Mmax •'max 



denotes the arrival time of the reflected PsPs head wave and the complex functions 
V := v{t) and 7 := 7(t) are defined by 



v{t) = -i 



y + h 1 t"^ xt 



y 7-— I— ^ lY* Zj iv* ^ 



(,\ ■ y±Ji /I t^ xt 



T \ T/+ r r 
V ^Ps 



for th < t < to and x < 0, 
for th < t <to and x > 0, 



and 



X y + h t^ 1 

7(i) = -i^tH W^ -o fort>to. 



Ps 
• V^^g is the velocity of the reflected PsS wave and satisfies: 



'^PssJ^'V^t) 



^PsS,y(^'y't) 



■ Ps 

IT 



^e 



F 



Ps 



vr 



^e 



K+{v{t))np,s{vit))—{t) 



dv 



■iv{t)npss{v{t))—{t) 



ift}^<t<to 

and \Qm (f(to))| > 



l^,r 



^pssA^^y^^") 



^PsS^yi^^y^^) 



F+ 
-^^^e 



vr 



F, 



Ps 



■K 



3f?e 



4(7(t))^Ps5(7W)^(i) 



d'y 



\i{t)npss{i{t))^{t) 



ift>to 



v^ Jx,y,t) =0. 



else 



Here to denotes the arrival time of the reflected PsS wave and tu denotes the arrival 
time of the reflected PsS head wave: 



h = y\ 



1 



V o ^ in 



' +/..M 



1 \x\ 

+ 



T/+2 y2 ' y 

V p "^ max ►max 



For tj^<t < tQ, the function v{t) is implicitly defined as the only root of 



1 



1/2 



1/2 



(?GC^^(g,t)=y — ^ + (7') +h 



V+^ 



+ > 



+ iqx — t, 



S / \ "Ps / 

such that Qm ( x "^^ ' I < 0. 

For t > tQ, the function 7(t) is defined as the only root of q G C ^^ J^{q,t) whose real 
part is positive. 
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I'pfPf is the velocity of the transmitted PfPf wave (the Pf transmitted wave generated 
by the Pf incident wave) and satisfies: 



^PfPf,xi^^y^i) 



^PfPf.yi^^y^t) 



-PiiFp. 



Pf 



^e 



vr 



du 
iv{t)TpfPf{v{t)) — {t) 



V7,F, 



IT 



^e 



KpAv{t))TpfPf{v{t)) — {t) 



dt 



iftf^<t<to 

and |9m (f(to))| > 



K. 



''pfpfj^^y^'t) 



V7,F. 



IIVP/ 

vr 



^e 



KiFpf 
^PfPf,y(.^^y^^) = ^e 



,jit)Tpfpf{j{t))^{t) 
i^-pf{l{t))TpfPf{-i{t))'^{t) 



dt 



ift>to 



i^PfPf{x,y,t) = 0. 



else 



Here to denotes the arrival time of the transmitted PfPf wave and tf^ denotes the 
arrival time of the transmitted PfPf head wave: 



1 1,1 1 \x 

th = ~y' I — ^ ~ "^ — ^^' / — ~ ~ -^^ — ^ 



V, 



Pf 



2 y2 



V p r max ^ max 



For tj^<t < to; the function v{t) is implicitly defined as the only root of 



1/2 



1/2 



qeC^ F{q, t) = -y\ ^-^ + g^ ) + /i ( ^^-^ + g^ ) ^ iqx - t, 



V, 



Pf 



V, 



Pf 



such that 9m ( x ^^ ) < 0. 



For t > to, the function j{t) is defined as the only root of q €z C ^^ J-{q,t) whose real 
part is positive. 

i^PfPg is the velocity of the transmitted PfPs wave and satisfies: 



''pfpsj^^y^*) 



^12-^P/ 



^e 



TT 



dt) 

yv{t)TpjPs{v{t)) — {t) 



— 12 P f 

^pfPs,y(^^y^''^) = '^^ 



dzj 
'^PsHi))'^PfPsi^it)) — {t) 



iftj^<t<to 

and 1 9m (t;(io))| > 



V„ 



^pfPs^xi^^y^t) 



v;oF, 



12VP/ 
TT 



3fie 



rv.F^, 



'^pfPs,yi^^y^'t) = — —^(^ 



TT 



vi{t)TpfPs{^{t))^{t) 
d'y 

^^Psi^{t))TpfPs{l{t))^{t) 



iftyto 



PfPs 
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^PfPs{x,y,t) = 0. 



else 
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Here to denotes the arrival time of the transmitted PfPs wave and tu denotes the arrival 
time of the transmitted PfPs head wave: 



1 1,1 1 \x\ 

tfi= -yj T - T7^ h/l, / — - T7^ h 






V pr '^ max "^ ir 



For tj^<t < tQ, the function v(t) is implicitly defined as the only root of 

( 1 V" ( 1 \"^ 



V, 



Ps 



V, 



Pf 



such that Qm ( x ^^^ ' J < 0. 

For t > to; the function ^{t) is defined as the only root of q & C ^^ J-{q,t) whose real 
part is positive. 

i^Pfg is the velocity of the transmitted PfS wave and satisfies: 



''pfS,xi^^y^*) 



'^pfsj^^y^'^) 



F. 



+ 

Pf 



vr 



^Pf 



vr 



3?e 



^e 



^-{v{t))Tpfs{vit)) — {t) 



dij 
vv{t)Tpfs{v{t)) — {t) 



iftf^<t<to 

and \Qm (f(to))| > 



Vr, 



^pfs,xi^^y^*) 



'^pfsj^^y^'^) 



F+ 
^Pf 



IT 



T\ 



Pf 



n 



^e 



^e 



^-s{i{t))rpfs{i{t))^it) 



i^{t)Tpfs{l{t))^{t) 



ift>to 



t^Pfsi^^y^t) = 0. 



else 



Here tg denotes the arrival time of the transmitted PfS wave and t^ denotes the arrival 
time of the transmitted PfS head wave: 



1 



^h = -y\ 



T/-2 y2 

V a max 



' +/. ' ' 



1 \x\ 

+ 



Vpp "^max ►max 



For tu<t<tQ, the function v{t) is implicitly defined as the only root of 



1/2 



q£ C^J'{q,t) 



V, 



+ 



+ h 



V, 



+ q +iqx -t, 



Pf 



such that Qm ( x "^^ ' I < 0. 

For t > tQ, the function 7(t) is defined as the only root of q G C ^^ J^{q,t) whose real 
part is positive. 
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^PsPf ^-^ ^^^ velocity of the transmitted PsPf wave and satisfies: 



'^PsPfJ^'y^i) 



PsPf,y 



-^ii^JRe 



TT 



x,y,t) 



7?- 77+ 



TT 



\v{t)TpsPf{v{t)) — {t) 



dv 



Kpf{v{t))TpsPf{v{t)) — {t) 



dt 



iftj^<t<to 

and 1 9m (t;(to))| > 



Kr 



^psPf,xix^y^i) 
^^psPfJ^^y^t) 






TT 



vr 



i7(t)rp,p/(7(t))^(t) 

d'y 
>^'pMt))TpsPf{l{t))-^{t) 



dt 



ift>to 



^PsPfi^iVi'^) = 0- 



else 



Here to denotes the arrival time of the transmitted PsPf wave and t^ denotes the arrival 
time of the transmitted PsPf head wave: 



th = -y^ 



1 



Vp^ ^^ax 



' ^h^ ' 



1 \x\ 

+ 



VL' ^^ax V, 



' Ps 



For ti^ < t < to, the function v{t) is implicitly defined as the only root of 



1 



1/2 



1 



1/2 



q e C ^ J'{q,t) = -y \ ^-ij + q^ ] + /i — n + « + iqx - t 



Vr 



_ 2 
Pf 



V, 



Ps 



such that 9m ( x "^^ ' I < 0. 

For t > tQ, the function ^{t) is defined as the only root of q £ C i-^ J^{q,t) whose real 
part is positive. 

• t^PsPs ^^ ^^^ velocity of the transmitted PsPs wave and satisfies: 



^PsPs,xi^^y^t) 



■^i^^Ke 



n 



TT 



dij 
ivit)TpsPs{vit)) — {t) 



dx) 

>^Ps{<t))'^PsPs{v{t)) — {t) 



ift}^<t<to 

and 1 9m {v(to))\ > 



► m 



^PsPsJ^^y^t) = -'^^^^^e 



TT 



^PsPs,y{^^y^t) = '^^^^e 



i^{t)TpsPs{l{t))'^{t) 
K-^{j{t))TpsPs{l{t))^{t) 



iftyto 



^PsPsi^^y^''^) =0- 



else 
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Here to denotes the arrival time of the transmitted PsPs wave and tu denotes the arrival 
time of the transmitted PsPs head wave: 



^h = -y\ 



T/- 2 y2 
V p^ '"^^^ 



' +/..U 



1 \x\ 

+ 



T/+2 y2 ' Y 

V p * max "^ max 



For t/j < t < to; the function v{t) is implicitly defined as the only root of 



1/2 



qeC^ F{q, t) 



such that 9m ( x ^^ ) < 0. 



_2 
Ps 



V, 



+ 



+ h 



V, 



2 V'' 
+ q +tqx -t, 



Ps 



For t > to, the function ^(t) is defined as the only root of q G C ^^ J-{q,t) whose real 
part is positive. 



• I'pgg is the velocity of the transmitted PsS wave and satisfies: 



^PssJx^y^^) 



^pss,y(^'y'0 



^^PssJ^^V't) 



'^Pss/^'V'*) 



F^s 


3f?e 


■K 


- 


^Ps 


^e 


vr 


L 


Fk 


^e 


vr 


_ 


•^Ps 


^e 


vr 





Kg{v{t))Tp,s{vit)) — {t) 



\v{t)rpss{v{t))'^{t) 



^shitWPss{l{t))^{t) 



ij[t)Tpss{i{t))^{t) 



■ifti^<t<to 

and |9m (f (to))| > 



K. 



ift>to 



^Pssi^^V^'t) =0- 



else 



Here to denotes the arrival time of the transmitted PsS wave and tu denotes the arrival 
time of the transmitted PsS head wave: 



1 



^h = -y\ 



Vf' v^ 



' ^h^ ' 



V, 



1 X 

+ 



+2 y2 ' Y 



max " max 



'S """^ V '^Ps 

For tj^<t < to, the function v{t) is implicitly defined as the only root of 



G 



^{q,t) = -y 



1 



1/2 



V. 



+ q' 



+ h 



1 



V, 



+ q +tqx -t, 



Ps 



such that 9m (x ^^^ I < 0. 

For t > to, the function ^{t) is defined as the only root of q £ C ^^ J-{q,t) whose real 
part is positive. 
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Remark 2.1. For the practical computations of the displacement, we won't have to explicitly 
compute the primitive of the velocities v, which would he rather tedious, since 

/ u{T)dTJ * f = iy* I f{T)dT 

Therefore, we'll only have to compute the primitive of the source function f . 

The proof of this theorem is similar to the one detailed in [8] for the computation of the 
analytical solution to the acoustic/poroelastic problem, therefore we won't detail it here. 

3 Numerical illustration 

To illustrate the use our results, we have compared our analytical solution to a numerical one 
obtained by C. Morency and J. Tromp [13]. We consider an two-layered poroelastic medium 
whose characteristic coefficients are 

• the solid density: pf = 2200 kg/m^ and p^ = 2650 kg/m^; 

• the fluid density: /ot = 950 kg/m^ and pT = 750 kg/m^ ; 

• the porosity: </>"'" = 0.4 and (j)~ = 0.2 ; 

• the tortuosity: a^ = 2 and a^ = 2; 

• the solid bulk modulus: K^ = 6.9 GPa and K^ = 37 GPa; 

• the fluid bulk modulus: Kf = 2 GPa and K^ = 1.7 GPa; 

• the frame bulk modulus: K^ = 6.7 GPa and K^^ = 2.2 GPa; 

• the frame shear modulus p'^ = 3 GPa and p^ = 4.4 GPa; 
so that the celerity of the waves in the poroelastic medium are: 

• for the fast P wave, Vpr = 2692 m/s and Vpr = 2535 m/s; 

• for the slow P wave, Vp^ = 1186 m/s and Vf^ = 744 m/s; 

• for the tp wave, Vg = 1409 m/s and Vg = 1415 m/s. 

The source is located in the acoustic layer, at 500 m from the interface. We used two types 
of sources in space: the first one is a bulk source such that fu = fw = —10^° and fp = 0; the 
second one is a pressure source such that fu = fw = and /p = 1. In each case we used a 
fifth derivative of a Gaussian of dominant frequency /o = 15 Hz: 



7r2 



fit) = 4 „2 
Jo 



l\ .vr2 / iV .vr4 / 1 



'''-tJ^'ts['-tJ -'ts[' h 



"^(*^i) 



for the source in time. We compute the solution at two receivers, the first one is in the upper 
layer, at 533 m from the interface; the second one is in the bottom layer, at 533 m from the 
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interface; both are located on a vertical line at 400 ni from the source (see Fig. [3]). We repre- 
sent the y component of the displacement from t = Otot = lsin Fig. |4] for the bulk source 
and in Fig. \5\ for the pressure source. The left pictures represents the solution at receiver 
1 while the right pictures represents the solution at receiver 2. On all the pictures the blue 
solid curve is the analytical solution and the red dashed curve is the numerical solution. 

All the pictures show a good agreement between the two solutions, which validates the nu- 
merical code. 







Receiver 1 


n+ 


Source 


' 


533 m 


400 m 
500 m 






■ 
















n- 






533 m 






• 

Receiver 2 



Figure 3: Configuration of the experiment 




PsPs 




Figure 4: The y component of the displacement at receiver 1 (left picture) and 2 (right 
picture) in the case of a bulk source. The blue solid curve is the analytical solution computed 
by the Cagniard-de Hoop method, the red dashed curve is the numerical solution. 



4 Conclusion 

We provided the complete solution (reflected and transmitted wave) of the propagation of 
wave in a two-layered 2D poroelastic medium and we used it to validate a numerical code. 
In a forthcoming paper we will use this solution as a basis to derive the solution in a three 
dimensional medium. 
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-Cagniard de Hoop 
- Numerical 



^PsPf 
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PfPf 




Figure 5: The y component of the displacement at receiver 1 (left picture) and 2 (right picture) 
in the case of a pressure source. The blue solid curve is the analytical solution computed by 
the Cagniard-de Hoop method, the red dashed curve is the numerical solution. 
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